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cn '. Abstract 
O 

' We study the perturbations on both "radial" and "angular" components of the quintessence with 

(N 

an internal abelian and non-abelian symmetry. The properties of the perturbation on the "radial" 

^3 : 

, component depend on the specific potential of the model and is similiar for both abelian and non- 

OO I abelian case. We show that the consine-type potential is very interesting for the 0{N) quintessence 

model and also give a critical condition of instability for the potential. While the properties of 
perturbations on "angular" components depend on whether the internal symmetry is abelian or 
I non-abelian, which we have discussed respectively. In the non-abelian case, the fluctuation of the 

00 

' "angular" component will increase rapidly with time while in the abelian case it will not. 

o 
o 
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1. Introduction 

Recent observations show that our Universe is fiat and 73 percent of its total energy 

density is resulted from "dark energy", which has a negative pressure and can accelerate the 
expansion of the Universe 0,0. using recent announcement by the WMAP satellitei, we 
still do not judge whether dark energy is due to an unchanging, uniform and inert "vacuum 
energy" (also known as a cosmological constant) or a dynamic cosmic field that changes with 
time and varies across space (known as quintessence). Quintessence can be considered as a 



spatially homogeneous scalar field that evolves with time and many mo dels [61,0,0, 9, 10, 11 1 
have been constructed so far. The possibility of quintessence as rolling tachyon has also 
been investigated [l^ . 

Boyle et al|l3| and Gu and Hwang|il4i] have discussed quintessence with complex scalar 
field. In a previous paper we have further generalized their ideas by replacing the 
complex scalar field with a A/^-plet scalar field which is spinning in a 0(A'^)-symmetric po- 
tential. When one of the angular components is fixed, this 0{N) quintessence model will 
reduce to the 0{N-1) quintessence model. If the N-2 angular components are fixed, 0{N) 
quintessence model can reduce to the complex scalar model mentioned above. If all angular 
components are fixed, 0{N) quintessence model will reduce to the quintessence. 

It is worth noting that this generalization does not hold its importance for the widely 
studied tracker-type potentials because the amplitude of the scalar field will increase steadily 
and make the angular contribution negligible. But for another type widely investigated 
potential, the cosine- type potential^, it will be very interesting. It has been pointed 
out that a natural introduction of quintessence is an ultralight axion with an almost massless 
quark|i3] and the cosine-type potential can be derived from such a model. Furthermore, This 
potential requires that the amplitude of the field should not goes to be very large and 
therefore, unlike the tracker-type potential, make the angular contribution significant. 

It would be very interesting to study the behaviours of the field when it is perturbed. 
We firstly investigate the fluctuation of the "radial" component and find that its stability 
depends on the specific potential of the model while has nothing to do with whether the 
the internal symmetry is abelian and non-abelian. We give a critical condition of the insta- 
bility for the "radial" perturbation, under which the "radial" fluctuation will grow rapidly 
with time. Fortunately, we find that this condition is not satisfed by most potentials for 
quintessence models. While the fluctuation of the "angular" components depend heavily 
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on whether the internal symmetry is abehan or non-abehan. In the abehan case (N — 2), 
the "angular" fluctuation will always damp, but in the non-abehan case{N > 2), there are 
possibilities for the " angular" perturbations to grow rapidly with time and become strongly 
space-dependent. It is worth noting that when we take the metric fluctuations(which are 
generally small compared with the fluctuations of the quintessence held) into account, the 
above conclusions are still held. In order to make the discussion more clear, in this paper 
we will restrict ourselves to the N — 2 and N — 3 cases which represent the abelian and 
nonabelian symmetry respectively. 

2. The Model 

We start from the flat Robertson- Walker metric 

ds"" ^ dt'' - a\t){dx^ + dy"" + dz'') (0.1) 
The Lagrangian density for the quintessence with 0(3) symmetry is 

= Ig'^id^^nidu^") - ^(1*1) (0-2) 
where is the component of the scalar held, a = 1,2,3. To make it possess a 0(3) 
symmetry, we write it in the following form 



= R{t) cos If i{t) 
<I>^ = R{t) smipi{t) cos (p2{t) 

= R{t) sin (/5i(t) sin (^2(^) (0-3) 

Therefore, we have |$"| = R and assume that the potential of the 0(3) quintessence depends 
only on R. It is clear that when we set the component (p2 to zero, the above 0(3) system 
will reduce to the 0(2) abehan case. 

The Einstein equations and equations of motion for the scalar flelds can be written as 

= )^ = . ^) + V(«)] (0.4) 
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(0.6) 



R 

Lpi + (3iJ + 2— )</3i - sin Lpi cos (pi(f2^ = (0.7) 
K 

R 

ip2 + {3H + 2-)ip2 + 2cotipiipiip2 = (0.8) 
R 

where H is Hubble parameter and the term^^p^ comes from the first integrals of the equations 
of motion for the angular components(detailed discussion seejl^). 

3. General Equations of Motion for "Radial" and "Angular" Perturbation 

In order to eliminate the ambiguity from gauge freedom, one has to identify gauge in- 
variant quantities or choose a given gauge and perform the calculations of perturbation in 
that gauge. In this paper, we will carry out our investigation in synchronous gauge just as 
Ratra and Peebles have done in Ref . . As we shall see, in synchronous gauge the evolution 
of the modes that we consider here depend on the gravitational field only through a gauge 
invariant quantity |, so we can consider a single perturbation mode conveniently in the 
following sections. 

The line element of perturbed metric of a spatially fiat FRW spacetime is taken as 

ds"^ = df - a'^{t){5ij - hij)dx'dx^ (0.9) 

where hij are the metric fiuctuations and | hij |<< 1. For simplicity, we write the first-order 
equations of perturbations for the N = 3 case. Decomposing the components of the field as 
follow 

R{t,x) = R{t) + 5R{t,x) (0.10) 



ipi{t,x) = ipi{t) + 5ipi{t,x) (0.11) 

ip2{t,^) = V2{t) + SMt,^) (0.12) 
We can obtain the equations of motion for the fiuctuations as 
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{5R) - - [ipi^ + sin" ipi]{5R) + 3-{5R) 

+V"{R){6R) - UiR - 2^iR{6^i) 

-Rsm{2ipi)<f2^{6^i) - 2Rsm'^{(pi)^2{S^2) = (0.13) 



(5^i) - 4v'(M) + 3-(M) + 2^('^^i) 
a a, JrC 

-(^2^cos(2(/?i)(5(/7i) - ^hifi - sin{2(pi)(p2{Sip2) 

+|v3i(5i?) - ^RM^R) = (0.14) 



{S'<f2) - ^ V'(5(^2) + 3-(5(^2) + '2^{S^2) - lhip2 

+2cot(v9i)[(/3i(5(^2) + '^2iSipi)] - 2cSC^{ipi)(^i(f2{S(fi) 

+^ip2{S'R) - ^R^SR) = (0.15) 



h + 2Hh^2R{SR) + 2R'^ipi{6ipi) + 2R^sm^{ipi)ip2{S^2) 

+ 2^(5i?) + R^sm{2^,)ip2\S^i) - V'{R){6R) (0.16) 



- ^ijj = Rdi{5R) + R^ipidi{5ipi) 

+ R^ sin\^^)ip2di{8ip2) (0.17) 



~{hij^kk + h^ij hij^j]^ hjf^ij^) 3Hhij 
a 

-Hh6ij - h'ij = 5ijV\R){5R) (0.18) 

where S^pi and 5ip2 are fluctuations of the " radial" and "angular" components re- 
spectively. The above equations of motion for fluctuations are the most general case for 
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quintessence with 0(3) internal symmetry. When we consider only the "radial" perturba- 
tion, we set ip2, 5lpi and to zero. If we deal with the perturbation on the "angular" 
component in the abelian case, what we need to do is to set another "angular" component 
ip2 together with its fluctuation 5lp2 to zero. In the subsequent sections, we will discuss them 
respectively. 

4. Perturbation on the "Radial" Components 

As we have pointed out in the introduction, the property of the perturbation on the 
"radial" component depend only on the potential of the model and is independent of whether 
the internal symmetry is abelian or non-abelian. In this section, we will show this in detail. 
As we all know that scalar fields with an internal symmetry are likely to produce Q balls 
or other non-topological solitons, which have been studied by many authors 16, Il8 |. 



In their work, they generally put the gravitational effects aside because the magnitude of 
gravitational fluctuations induced by the fluctuations of the scalar fields are far smaller 
than the self-interaction of the scalar fields. Firstly, we, following the previous study in this 
field|l6, 17, l^. 

also investigate the perturbation without considering the metric fluctuation. 



After this, we take the gravitational effects(metric fluctuations)into account and carry out 
a similar study on the fluctuations. The equations of motion for fluctuations of the "radial" 
component are 

+V"{R){5R) = (0.19) 



{6R) + 3H{6R) - i^V'{6R) + -^,m) 



which is obtained by setting 5ipi, bip2 and h in Eq. ()0.13p to zero. If we choose for the 
fluctuation the following form: 

bR{t, x) = <5i?o exp[a(t) + ^kx] (0.20) 
then for nontrivial bR^ , we have 

Following the authors in Ref . (ifil. Il7l |. we assume that a(t) is a slow- varying function, i.e. 
ait) « d? and a ~ Const. In the following sections, we shall always hold this assumption. 
Therefore, neglecting the a in the above equation, we have 



1 

"=2 



(0.22) 



If a is real and positive, the fluctuation will grow rapidly with time. Therefore the 
instability band for this fluctuation is 

From Eq. ()0.23p . one can find that the instability band depends on the specific potential of 
the quintessence. 

When considering the metric fluctuations, one can obtain the following equations of 
motion for the fluctuations 



+V"{R){6R) -Hr = (0.24) 



Q2 

h + 2Hh = 2R{6R) + 2^^(5^) - V'{R){6R) (0.25) 



h^i - hijj = Rdi{SR) (0.26) 



-Hh5ij - h,j = 6ijV' (R) {6R) (0.27) 
If we choose = 0, i.e. in the case of = 1 the Eqs. ()0.24|) - ()0.27p will reduce to Ratra- 



Peebles's results in the absence of baryonic term ^. Clearly, since the equations of motion 
for the metric and scalar fluctuations(Eq.( )0.25j) and Eq. ()().24|l ) are linear equations, the 
fluctuations could be taken as the following form 

6R{t, x) = 6Ro exp[a(t) + zkx] (0.28) 
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h{t,x) = hoexp[a{t) + zkx] (0.29) 

Since there are no V^/i term in Eg. ()().25j) . h can not oscillate rapidly and in fact, h will 
not be able to react, in lowest order, to the rapidly oscillating source terms in 
Then for nontrivial SRq and h^, we have 



a + 3H — --rr^R I a — 



2 6R 







a 



^ a2 ^ dR^ 







So 



a 



(0.30) 



(0.31) 



±W(3^-i^^)^-4(^ 



+ 



2 6Ro"' ^^a? a^R'^ dR^ 
From Eq. ()0.3ip . it is clear that the instability band is the same as that when we did not 
consider the metric fluctuations. This is reasonable in physics because the metric fluctuations 
induced by the fluctuations of the field are far smaller than the fluctuations of the field and 
its back-reaction to the field would be even smaller and thus negligible. So, when considering 
the metric fluctuations, there should not be a substantial change on the properties of the 
fluctuation of the quintessence field. 

5. Perturbation on the "Angular" Components in the Abelian Case 

In this section, we investigate the perturbation on "angular" component in the abelian 
case, that is, the case in which the quintessence fields possess a 0(2) internal symmetry. 

By setting 6R, h, 6ip2 and in Eq. ()0.14p to zero, we can obtain the equation of motion 
for the fluctuation of " angular" component up to the first order as: 



S'^i + {3H + 2^)6^1 - \v^5ipi = 
R 

where dipi is the fluctuation of the "angular" component. If we choose 

Sipi{t,x.) = 6(/)ioexp[a{t) + ikx] 
then for nontrivial Sipio, from Eq. ()0.32|) we have 



(0.32) 



(0.33) 



+ {3H + 2-)a + — = 
R 



(0.34) 
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It is clear that 



a 




(0.35) 



From Eg. ()().35j) . a is always negative and thus the fluctuation will damp quickly. That is, 
the angular perturbation will not produce a significant "angular" inhomogeneity and the 
global symmetry will not likely to become space-dependent. 

In the following, we will introduce the metric fiuctuations into our analysis. By a similar 
procedure as that in last section, we can obtain the equations of motion for the fiuctuations 
of metric and "angular" component up to the first order as following: 



^'h + Hh = R'^ipidipi 



(0.36) 
(0.37) 



2 {hij,kk ~\~ h^ij hikjk hji^^i}^ 3Hh 



-Hh5in — h 



(0.38) 



(5^1 + {3H + 2^)5^1 - 4 V'^v^i - Ihipi = (0.39) 

Since Eq. ()0.36|) and Eq. ()0.39|) are linear equations, the fiuctuations could be chosen in the 
following form 



5ipi{t,x) = SLpioexp[a{t) + ^kx] (0.40) 

h{t,x.) = hoexp[a{t) + zkx] (0.41) 
Then for nontrivial 6lpio, we have (from Eq. ()0.39p ) 

+ {3H + 2- - --^ip,)a + — = (0.42) 



R 2 6(^10 d 
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Therefore, one can obtain 



a 



R 1 ho 



(0.43) 



From Eq. (j0.43p . it is clear that a will always be negative even if the metric fluctuations are 
considered. 

6. Perturbation on the "Angular" Components in the Non-abelian Case 

In this section, we generalize the discussions in last section to the case in which the 
internal symmetry is non-abelian, that is the symmetry group is 0(3). We restrict ourselves 
to the case that only one "angular" component is perturbed. This will not lose its generality 
but greatly facilitate the discussion because we can always choose a coordinate system in 
which the perturbation appears in one angular direction. By setting 6R, 6(f2 and h in 
Eq. (j0.14j) to zero, we can obtain the equation of motion for the angular fluctuation up to 
the first order as following: 



S'^, + {3H + 2^)6^1 - 



R' 



V 5(pi — cos2(/?i(/32 Sifi = 



(0.44) 



where (fi and ip2 are the homogeneous parts of the "angular" components. If we choose for 
6lpi the same form as in Eq. (j0.40p . then for nontrivial SifiQ, from Eq. (|0.44j) we have 



R 
R 



cos 2ipi(^2^ = 



It is clear that 



(0.45) 



a 



(3 J/ + 2^) ± W (3ff + 2^)2 - 4 ("^ - cos 2^iv52' 



From Eq. ()0.46p . one can find that a could be positive if 



— cos 2(/7i(/32 < 



(0.46) 



(0.47) 



That is, under the above condition fEq . ()0 . 47|) ) the "angular" inhomogeneity might grow 
rapidly with time and make the symmetry group become space-dependent. 
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When taking into account the metric fluctuation and following a similar process as that in 
section 5, we can obtain the equations of motion for the fluctuations of "angular" component 
as follow 



R 



1 . 



— cos2(/?i(/32 6(fi hifi = 



(0.48) 



The equations of motion for metric fluctuations are 



(0.49) 



(0.50) 



ik 



3Hh 



-Hh5in — h 



(0.51) 



Similarly, we choose the fluctuations to be the form of Eq. ()0.40|) and Eq. ()0.4ip and for 
nontrivial SipiQ, from Eq. ()0.48|) we have 



+ {3H + 2^- l^ipi)a - cos 2ipiip2^ + ^ = 
R 2 Oipio 



and therefore 



a 



R 1 ho 



R 1 hn f k'^ r. 



(0.52) 



(0.53) 



From Eq. (j0.53|) . one can easily identify that the condition for positive a is the same as that 
in Eq. (|0.47p . This shows that in the 0(3) case, the instability condition for the "angular" 
fluctuations won't be changed when considering the metric fluctuation, just as that we have 
proved in the 0(2) case. 
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7. Conclusion And Discussion 

In this paper, we investigate tlie perturbation on both "radial" and "angular" components 
of the quintessence fields with an internal abelian and non-abelian symmetry. We find 
that the fiuctuation of the "radial" component depends on the specific potential of the 
quintessence model. Under certain condition, this fiuctuation could grow rapidly and thus 
make the "radial" component space-dependent. For the "angular" perturbation, we find 
that the properties of the fiuctuations depend on whether the internal symmetry group is 
abelian or non-abelian. In the abelian case, the fiuctuation of the "angular" component will 
damp rapidly with time and therefore the symmetry group will not become space-dependent. 
In the case that the internal symmetry is non-abelian, we find that under certain condition, 
the "angular" inhomogeneities might increase rapidly with time and make the symmetry 
group space-dependent, or local. Here we briefiy interpret the physical meaning of this 
instability: 5ipi represents a small internal angular inhomogenity; This inhomogenity will 
increase rapidly under the conditionfEq. ()().47|l ). Note that this angular inhomogenity is not 
that of space-time. It is surely interesting to study the quintessence with a "local" internal 
symmetry, which we will investigate in a preparing work. When taking into account the 
metric fiuctuations induced by the fiuctuations of the quintessence field and assuming that 
the back-reaction of the metric fiuctuations on the quintessence field are negligible as shown 
in Ref. 0], the above conclusion still hold true. 

It is worth noting that we choose the non-abelian symmetry group as 0(3) in this paper. 
It is not difficult to generalize the 0(3) case to 0{N) case and one may find that in the 
0(A'^) case, the above conclusion for the non-abelian symmetry group still hold true even if 
the specific condition under which the inhomogeneity increase will change. 
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